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ABSTRACT 

We propose a simple theoretical model for star formation in which the local star formation rate in 
a galaxy is determined by three factors. First, the interplay between the interstellar radiation field 
and molecular self-shielding determines what fraction of the gas is in molecular form and thus eligible 
to form stars. Second, internal feedback determines the properties of the molecular clouds that form, 
which are nearly independent of galaxy properties until the galactic ISM pressure becomes comparable 
to the internal CMC pressure. Above this limit, galactic ISM pressure determines molecular gas 
properties. Third, the turbulence driven by feedback processes in GMCs makes star formation slow, 
allowing a small fraction of the gas to be converted to stars per free-fall time within the molecular 
clouds. We combine analytic estimates for each of these steps to formulate a single star formation law, 
and show that the predicted correlation between star formation rate, metallicity, and surface densities 
of atomic, molecular, and total gas agree well with observations. 

Subject headings: galaxies: ISM — ISM: clouds — ISM: molecules — stars: formation 



1. INTRODUCTION 

The last decade has seen a revolution in our under- 
standing of star formation in galaxies, driven by the 
advent of spatially resolved multi- wavelength surveys. 
Prior to this work, our observational constraints on the 
star formation process were largely limited to low res- 
olution surveys that characterized entire galaxies using 
only a handful of observable quantities, e.g. the mean 
surface density of star formation averaged over a whole 
disk. While these surveys yie lded a number of intriguing 
results - most famously the iKennicutH ()1998[ ) star for- 
mation law - they left unanswered many basic questions 
about the physics of star formation. For example, they 
could not clearly determine whether star formation cor- 
relates more stron gly with the molecular or total gas con - 
tent of a galaxy (iKennicuttI [19981 IWong fc Blitg l2002f ) . 
or whether star formation is regulated primarily by lo- 
cal processes within indivi dual star-f orming c louds (e.g. 
Krumholz &: McKed [20051 hereafter iKMOSi : IShu et al.l 



rable resolution and accuracy (iGil de Paz et al.l[200l . 
Observations of H i emission from the V LA as part of 
The H I Nearby Galaxy Survey (THINGS; IWalter et all 
I2008D and of CO emission by the BIMA Surve y of 
Nearby Galaxies (BIMA SONG: Heifer et all [20031) and 
the HERA CO-Li ne Extragalactic Survey (HERACLES; 
iLerov et al.]l2009l ) using the 30 m IRAM telescope pro- 
vide maps of the gas content of galaxies at compara- 
ble resolutions. Combining these data sets leads to two 
particularly important conclusions that our theoretical 
models of star formation must incorporate and explain. 

The first observational result is that that star for- 
mation is a direct product of the molecular gas in a 
galaxy, not of all the gas. Across a wide range of 
galactic environments the star formation rate corre- 
lates well with molecular gas, and poorly or not at all 
with the atomi c gas measured on sub -kpc spatial scales 
JWonK & Blitz 2002; 'K ennicutt et al.ll2007l : fLerov et all 
[200 8: Bigiel ct al, 2008i, hereafter 



20071 ) or by galactic-scale processes such as s piral shocks , 
supernovae, or cloud -cloud interact ions (e.g. IWvsill986t 
[Silk 1997; .Tan 200(^ ITJeTaDlmi^h . 

Now, however, emission maps at 24 /im from the 
Spitzer Infrared Nearby Galaxy Survey (SINGS) pro- 
vide us with accurate estimates of the rate of dust- 
enshrouded star formation at resolutions of better than 
a kpc in nearby galactic disks (Kennicutt ct al. 2003, 
l2007t ICalzetti et al.ll2007l) , while ultraviolet observations 
from the GALEX Nearby Galaxy Survey (NGS) reveal 
the rates of non-obscured star formation with compa- 
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The correla- 
tion between star formation rate and molecular gas sur- 
face density is super-linear in samples that extend to 
starbursts with gas surface densities ~ 10'^ Mq pc^^ 
(jKennicuttI 119981 ) . but is nearly linear for galaxies with 
molecular su rface densities from 5 — 100 Mq pc~^ (BOl; 
however, see iKennicutt et al.l l2007l for a different inter- 
pretation). Observations of the low density outskirts of 
galacti c disks hint that the linearity may break dow n 
there (jCardan et al.l 120071: iFumagaUi fc Gavazzil [2001) . 
but it is unclear if this indicates a change in the star for- 
mation process or a change in the CO to H2 conversion 
factor. 

The second observational result is that giant molecular 
clouds (GMCs) have remarkably similar properties in all 
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nearby galaxies. Across the Local Group GMCs appear 
to have the same surface density, roughly 85 Mp^ pc~^, 
and to obey the sam e lincwi dth-size relation (|Blitz et al.l 
I2007t iBolatto et al . 2008: Hever et alll2009f) . Together 
these two observations imply that all observed molecular 
clouds are not far from virial balance between gravity 
and internal turbulence. There is of course considerable 
uncertainty in the GMC surface density, arising mostly 
from the CO-to-H2 conversion factor, but it is striking 
that there is no clear evidence of a systematic trend in 
GMC properties across a sample of galaxies ranging from 
H i-dominated dwarfs to molecule-rich giant spirals. This 
seems to be a critical clue to the physics of molecular 
clouds. 

Efforts to incorporate these observational results and 
their implications into numerical simulat i ons a re al- 
ready underway. [Robertson fc Kravtsovl (|2008[ ) and 
iGnedin et al.l (|2009f ) present calculations of galaxy evo- 
lution that include models for the chemistry of hydro- 
gen molecule formation and destruction, and that re- 
strict star formation to occur only in molecular gas. 
iTasker fc Ta3 (l2009l) use high resolution simulations of 
galactic disks to study the origin of the observed proper- 
ties of GMCs. Our goal in this paper is to complement 
and extend this work by developing a simple analytic 
model for the star formation law that is based on the 
insights provided by the new observational surveys, and 
that incorporates the theoretical understanding that has 
developed around them. 

Our basic approach is to break the problem into three 
tractable pieces, each of which can be treated using 
models and observations already available in the liter- 
ature. The first is the problem of calculating what 
fraction of the gas in a given portion of a galaxy 
will be in the molecular phase and thus eligible to 
form stars. The problem of determining the molecular 
fraction in the ISM has been treated extensively (e.g. 
van Dishoeck fc Blacklll986t ISternberej|1988t lElmeereenl 
1993t iDraine fc Bertoldilll996t iBrowning et al.l 120031): in 
our approach we will adopt the model of lKrumholz et al.l 
(j2GG8, 20Q3, hereafter KMT08 and KMT09), which 
agrees very well with the observations o f molecular 
fractions reported by iBlitz fc Rosolowskvl (j2004i ) and 
iLeroy et"al1 (|2008l ). and has the advantage that it does 
not depend on unknown and generally unmeasurable 
quantities such as the intensity of the ultraviolet radi- 
ation field inside a galaxy. The second problem is to 
estimate the characteristic properties of the GMCs in a 
galaxy, which can be done using a combination of their 
observed properties simple arguments ba sed on virial bal- 
ance and th e galactic Jeans mass (e.g. iKim fc Ostriked 
120011: IKM051. The third problem is to estimate the rate 
at which molecular clouds of known properties trans- 
form themselves into stars. This rate is known from 
observations to be ^ 1% o f the mass per free-fall t ime 
(jZuckerman fc Evans! [l97l iKrumholz fc TanI 120071 ). a 
value that can be understood theoretically as a result 
of regulation of star formation by supersonic tu rbulence 
(iKravtsovl 12001 IKM05I : iWada fc NormanI 120071) . In § H 
we explain how these three components can be combined 
to produce a star formation law, and then in §[3] we com- 
pare the results to observations. 

2. THE STAR FORMATION LAW 



Given that star formation occurs in molecular gas, we 
formulate our theoretical law for the local star formation 
rate (SFR) surface density S, in a galaxy as a product 
of three factors: 

i;, = ^g/H2 ^ ■ (1) 

Here Eg is the total gas surface density at some point 
in the galaxy. In practice this will always be an av- 
erage over some size scale, determined by the resolu- 
tion of observations (or simulations). This determines 
the total available "raw material" for star formation. 
The factor is the fraction of this mass in molecular 
form; atomic gas does not participate in star formation. 
The molecular component of the gas is organized into 
clouds which have some mean volume density /9h2 j s-nd 
= [37r/(32G/9H2)]^^^ is the free-fall time at this mean 
density. The quantity SFRg is the dimensionless star 
formation rate; it is the fraction of the gas transformed 
into stars per free-fall time. Alternately, one may think 
of it as the star formation efficiency over one free-fall 
time (as opposed to the total star formation efficiency, 
which might mean the fraction of gas transformed into 
stars over some other timescale, such as the galactic ro- 
tation time or the lifetime of an individual GMC). The 
third factor, SFRff/tff, is simply the SFR per free-fall 
time divided by the free-fall time, which is the inverse 
of the time required to convert all of the gas into stars. 
To make a model for the star formation law, we must 
estimate /hj, iff, and SFRff in terms of the observable 
quantities for a galaxy. 

2.1. The Molecular Fraction 

The molecular mass fraction is determined by 
the balance between dissociation of molecules by the 
far ultraviolet (FUV) interstellar radiation field in the 
Lyman- Werner bands and formatio n of molecul es on 
the surfaces of dust grains. IKMT08I and IKMT09I show 
that to good approximation within a single atomic- 
molecular complex is a function of the gas surface density 
of the complex Scomp and the mctallicity Z . We will not 
repeat the full derivation of this result here, but a sum- 
mary of the calculation that produces it is that one first 
solves the idealized problem of finding where the tran- 
sition between the atomic envelope and the molecular 
interior occurs within a uniform sphere of hydrogen gas 
and dust embedded in an isotropic dissociating radiation 
field. This analysis shows that the fraction of the com- 
plex in molecular form depends on the dust optical depth 
of the complex ScompCd and on the dimensionless ratio 
X oc crdGo/(ncNM^), where cfd is the dust cross section 
per hydrogen nucleus. Go is the intensity of the dissoci- 
ating radiation field, ncNM is the number density of gas 
in the cold atomic medium that surrounds the molecu- 
lar part of the cloud, and TZ is the rate coefficient for 
II2 formation on the surfaces of dust grains. Since Ud 
and TZ are both, to first order, simply measures of the 
total amount of dust in a galaxy, their ratio should not 
vary widely between galaxies. Similarly, in a galaxy with 
a two-phase atomic medium the cold atomic gas density 
'T'CNM is determined by thermal pressure balance between 
the two phases, which in turn depends on the balance be- 
tween heating by FUV photons and atomic line cooling 
in the atomic gas. Analysis of these processes implies 
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that the ratio Gq/jt-cnm is a weak function of metal- 
licity and is otherw ise independent of galaxy properties 
(iWolfire et al.l[200l . Thus X varies httle between galax- 
ies, and this result enables us to write the molecular frac- 
tion for a given atomic-molecular complex as a function 
solely of its gas surface density Scomp and its metallicity 
Z: 



/H2(^coinp, Z') ?» 1 



4 1 +(5 



-5 



-1/5 



(2) 



where s = ln(l + 0.6x)/(0.04Eeomp,o^'), X - 0.77(1 + 

3 2^/0.365)^ J ^ 0.0712(0.1s-l + 0.675)^2.8^ Ecomp,0 = 

Scomp/(lAf0 pc"'^), and Z' is the metallicity normal- 
ized to the solar value. Note that this a pproxima tion 
is slightly different with the one given in IKMTOQI : the 
two agree to within a few percent for clouds that are 
substantially molecular, but this one is more accurate at 
small molecular fractions (McKee, Krumholz, & Tumlin- 
son, 2009, in preparation). 

Here Scomp is the surface density of a ~ 100 pc-sized 
atomic-molecular complex. However, extragalactic ob- 
servations generally measure a gas surface density Sg 
that is averaged over a much larger scale. Since in- 
creases super-linearly with Ecomp, clumping of the gas on 
scales below the observational resolution would lead us 
to underpredict if we were simply to use the large- 
scale-averaged value of Eg in place of Scomp equation Q . 
Since we wish to propose a model that is applicable to 
data and simulations at a range of resolutions, it is con- 
venient to approximately correct for this effect by letting 
Ecomp = cEg , where c > 1 is a clumping factor and c — > 1 
as the resolution approaches ^ 100 pc. 

As a final cave at, it is im portant to poin t that our cal- 
culati on of in lKMTOgI assumes that the IWolfire et al.l 
(|2003D semi-analytic model for the atomic ISM is applica- 
ble, and the model begins to break down a t metallicities 
bel ow roughly 5% of solar (see Figure 13 of lWolfire"eran 
[200 3il because dust grains and polycyclic aromatic hy- 
drocarbons begin to be neutral rather than positively 
charged, as the model assumes. Turbu lent heating of the 
cold H I phas e (|Pan fc PadoanI [20091 ). which is not in- 
cluded in the IWolfire et al.l models, is also likely to be 
important at low metallicity. Thus, although our gen- 
eral method of calculating molecular fractions will apply 
even at low metallicities, the relationship between ncNM 
and Go which is used to derive equation ([2]) is not valid 
at metallicities Z' < 0.05. 

2.2. Giant Molecular Cloud Properties 

Next we must compute tff and SFRff, which will de- 
pend on the properties of the star-forming GMCs in 
a galaxy. Before proceeding with such a calculation, 
we note that observations of Local Group galaxies indi- 
cate that in galaxies ranging from metal-poor dwarfs to 
molecule-rich spirals, the molecular cloud surface den- 
sity Eel ~ 85 Mq pc~2 and the molecular cloud virial 

ratio ttyir ~ 2 independent of gal actic environment 

(jBlitz et al.ll2007l:lBolatto et al.ll2008l ).^ This invariance 

^ Note that this value of Sj-j is lower jthan the 170 Mq pc~^ 
for Galactic GMCs found by Solomon et al. 1983i but is consis- 
tent with the low er value determined by the more recent survey of 
IHever et aL|[2009l . 



is reasonably easy to understand on theoretical grounds. 
The virial theorem implies that the mean pressure 
within the cloud is Pd/fce = 0.7 x 10^ OvirE^,^ K cm"^ 
(.KMOa). where E', = E H/(85Mfn pc'^). In comparison, 
iBoulares fc Coxl (jl990D find that the mean kinetic pres- 
sure^ in the ISM of a Milky Way-hke galaxy is 1.4 x 10"* K 
cm~^, an order of magnitude lower, although the pres- 
sure may be higher than average in spiral arms where 
GMCs form. The pressure is almost certainly lower 
in low surface-density dwarfs. The mismatch between 
the pressures in GMCs and the pressures in their en- 
vironments indicates that external pressure is at most 
marginally important in determining the properties of 
molecular clouds, and that GMCs must instead be in- 
ternally regulated; in effect, a CMC forgets about its 
galactic environment. Moreover, such a picture provides 
a quantitative explanation for the observed values of the 
CMC surface density and virial ratio. The dominant 
mechanism of i nternal regulation is H 11 region feedback 
(|Matzneill200l l. a process whos e efficiency depends on 
the column density of the cloud. iKrumholz et al.l ()2006l ) 
find that H 11 regions stabilize GMCs at column densities 
of roughly 100 Mq pc"^ and virial ratios a^ir ~ 1 — 2, 
consistent within the uncertainties with the observed val- 
ues. 

The constant surface densities and virial ratios of 
GMCs provide a natural way to estimate ts and SFRg. 
Consider a CMC of mass M, surface density Ed, and 
virial ratio avir- The volume density in this cloud is 
(37ri/V4)E|^/^M-i/2, so the free-fah time is 



iff = 8E^7-''/X^' Myr (3) 

where Me = M/10^Mq. Similarly, iKMOSi show that 
the star formation rate per free-fall time in a turbulent 
medium is approximately 

, -0.68 A ^-0.32 



SFRff 



(4) 



where is the 1-D Mach number of the turbulence 
and ecorc is the fraction of the mass in a gravitationally- 
bound prestellar core that is incorporated into a star 
rather than being ejected by protostellar outflows. iKMOSi 
adopt Ecorc = 0. 5 based on analytic mo dels showing 
ecorc « 0.25-0.75 (jMatzner fc McKeell2000t) ; more recent 
work suggests the true value is ecorc « 0.3 (jAlves et al.l 
I2OOI . so we adopt ecorc = 0.3. 

The virial ratio is related to the one-dimensional 
velocity dispersion cr in a CMC by avir = 
57r-i/2(AfEci)-i/VVG (iBertoldi fc McKed [l99l . 

so (7 = 3.7 al(^I^'^/^ Mg^'^ km s~-^.^ For a molecular 
cloud temperature of 10 K the corresponding Mach 



number is A4 = 

SFRff « 0.017a 



20ati:KrM.' ,so 

0.84'r^/-0.08 7i.r-0.08 



(5) 



This is consistent with the observed value SFRff « 0.01 
(jKrumholz fc Tanll200l . 

^ We consider only turbulent and thermal pressure because the 
mean galactic magnetic field and cosmic rays pervade GMCs and 
the intercloud medium equally, and therefore provide neither sup- 
port nor confining pressure. 

^ For our fiducial Sj.j = 1 and Ovir = 2, this agrees with the 
observed linewidth-size relation a ^ 0A4tg\l(R/pc)°-^°'^°-^° km 
s-l lIBolatto et al.ll2008l) to within the error bars. 
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Combining equations Q and (O, and adopting a fidu- 



cial value of ttv 



gives 



SFRfi 



-0.33 



(6) 



iff 0.8 Gyr ' 

The invariance of molecular cloud properties that we 
observe in nearby galaxies must break down in galaxies 
with sufficiently high surface densities, where the exter- 
nal pressure is no longer negligible compared to a CMC's 
internal pressure. Since pressure varies as P oc Y? for 
both molecular clouds and galactic disks (.KM05 ). the 
galactic environment will become significant in determin- 
ing molecular cloud properties once the galactic surface 
density averaged over large scales becomes comparable to 
the surface density of an individual CMC. In this case the 
CMC surface density must increase in order to maintain 
pressure balance with the rest of the galaxy's ISM, which 
simply requi res that Sci ~ for > 85 Mq pc~^. 
(Alternately, iKomugi et al.ll2006l suggest that a change in 
CMC properties might be expected when Eg ^ 10^ — 10"^ 
Mq pc~^ because at such high surface densities collisions 
between GMCs become common.) Observations are con- 
sistent with this hypothesis: in the central kpc of M64, 
where the galactic surface density runs from ~ 50 — 1000 
Mq pc~^, the CMC surface density is not constant, and 
instead rises with galactic pressure. Averaged over the 
entire galaxy the mean CMC su rface density is 250 Mq 
pc-2 (|Rosolowskv fc Blitd[2005l) . Thus, the free-fall time 

in CMCs in high surface density galaxies varies as Eg 
If we adopt a column density of Ed = 85 Mq pc~^ for all 
CMCs in normal surface density galaxies and Ed = Eg 
at higher galactic surface densities, then we have 



SFRff 

iff 



0.8 Cyr 



1, 



85 Mq pc- 



0.67 



(7) 



Equation ([7]) gives an estimate for SFRff /iff in a molec- 
ular cloud of a known mass. To complete the calculation, 
we must estimate the c haracte ristic molecular cloud mass 
in a galaxy. We follow 'KM05' in estimating that this will 
be determined by the Jeans mass in the galaxy, which is 

rr! TT^G^E^Q* 

(8) 



M 



G2E, 



4rj4 



where (jg is the gas velocity dispersion, Q is the Toomre 
Q of the galactic disk, and Vl is the angular velocity 
of its rotation. If we can directly measure Eg, SI, and 
Q, or Eg and dg, for a galaxy, then we can solve for 
M directly and substitute into equation d?]) to obtain a 
characteristic value of SFRff /iff for that galaxy. How- 
ever, often one or more of the quantities are unknown, 
and even when they are known it is useful to have a 
rough estimate in terms of a single quantity such as Eg 
rather than three quantities Eg, fi, and Q. Since Mq 
enters the star formation rate only to the 0.33 power, 
any errors we make in this approximation are unlikely to 
have strong effects. We therefore follow KM0 5:' in assum- 
ing that all galaxies will be marginally Toomre stable, 
Q « 1, and noting that there is broad statistical corre- 
lation r^/Myr"^ « O.O54(Eg/85M0 pc-2)0-49. If we use 
this correlation in IHl) then we obtain 



Mfi « 37 



85 Mq pc- 



1.0 



(9) 



Finally, it is worth noting here that our estimate of the 
molecular cloud volume density, which de pends on Ed 
and Mg, is somewhat different than that of lKM05l They 
assumed that CMC surface densities were set largely by 
external pressure in a galaxy, and computed the den- 
sity based on this assumption. As discussed above, more 
recent observational and theoretical work suggests that 
instead CMC densities are primarily set by internal feed- 
back processes and do not vary significantly with galac- 
tic conditions, at least in Milky Way-like galaxies. Our 
model in this paper takes this result into account. 

2.3. The Full Star Formation Law 

We have now derived the major components of our star 
formation law (equation [1]) . The molecular fraction 
depends only on gas surface density Eg, metallicity Z' , 
and the clumping of the gas c on scales unresolved in a 
given observation or simulation (equation^). It increases 
with Eg, becoming fully molecular at ~ 10/cZ' Mq pc^^. 
We have also derived an analytic relation for the inverse 
star formation timescale SFRff /iff in two regimes. Where 
internal CMC pressure far exceeds the ambient ISM gas 
pressure and CMCs "forget" their environment - as typ- 
ically occurs in nearby galaxies with Eg < 85 Mq pc^^ 
- this timescale does not depend on Eg except indirectly 
through the molecular cloud mass (equation [9|) . Above 
Eg = 85 Mq pc~2, ambient pressure becomes compara- 
ble to the CMC internal pressure and the star formation 
timescale depends on Eg (equation [7]) . In neither case 
docs the timescale depend on either the metallicity or 
the clumping, so the star formation rate in molecular 
gas does not depend on either of these quantities. Only 
the star formation rate in total gas does. 

We are now ready to combine these pieces into our 
single star formation law: 



(10) 



3. COMPARISON TO OBSERVATIONS 

We compare our proposed star formation law, equa- 
tion (fTO|) . to the observed relationship between star for- 
mation, atomic gas, and molecular gas in Figures [T] 
and [21 The majority of the observations come from 
the THINCS sample. The full sample covers metallic- 
ities fro m logZ' = -1.22 to 0.49 (Walter ct al. 200i; 
IKMTO90 , but only four of the thirty- four galaxies have 
metallicities below \ogZ' = —1.0, and these are all 
dwarfs with such low star formation rates that they con- 
tribute negligibly to the total star formation rate in the 
sample. Moreover, the molecular gas masses for these 
systems are likely to be extremely uncertain (see below). 
Thus we adopt \ogZ' = —1.0 to 0.5 as a realistic range 
of metallicities in the data. 

The THINCS sample is observed at a resolution of 
~ 750 pc, much larger than a single atomic-molecular 
complex, so we expect c > 1. The true value of c 
cannot be determined directly in external galaxies with- 
out higher resolution observations. A lower limit comes 
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Fig. 1. — Star formation rate surface density E, as a function of H I {panel a), H2 {panel 6), and total gas {panel c) surface densities 
Sjji, Shji and Eg. Lines show our theoretical model predictions for values of clumping fact or tim es metallicity of logcZ' = —0.3,0.2,0.7, 
and 1.2, as indicated. Contours show observations from THINGS, and are constructed as in lBOSt we break the plane of the plot into bins 
0.05 dex wide in each direction and count the number of independent data points in each bin. The contours represent, from lightest to 
darkest, 1, 2, 5, and 10 data points. The dashed vertical lines in the Ejjj and Eg plots indicate the THINGS CO sensitivity limit of 4.5 
Mq pc~^. Note that our plots are shifted by a factor of 1.36 relative to those of IB08I because we include the mass of helium in Ehi, Ejjj, 
and Eg. 



from the fact that the observations mix together spiral 
arm and inter-arm regions, and the arm-interarm density 
contrast is ^ 2 — 4 in g alaxies observed at higher res- 
olution (jNakanishi fc So fue 2003; Schuster ct al. 2007). 
The complexes themselves represent density peaks on 
top of the already-enhanced density within the arm, and 
in fully molecular regions clouds are observed to have 
surface densities higher than the mean by a factor of 
2 (Rosolowsky & Blitz 2005). We therefore adopt 
c « 5, and thus we expect the data to be characterized 
by \ogcZ' « —0.3 to 1.2, with the four low-metallicity 
dwarfs lying at somewhat lower log cZ' . 

Our simple model recovers a number of salient features 
in the observations. Figure [T^ shows that we recover the 
observational result that the H i surface density reaches a 
maximum value, which is ~ 10 Mq pc~^ at Solar metal- 
licity, and that the star formation rate does not correlate 
well with Shi in resolved observations of galaxies. The 
RMS noise in the star formation rate surface densit y in 
the survey is - lO"'^ Mq kpc"^ (" Bigiel et al.]|2008D . so 
the apparent flattening of the contours below this value 
is an observational artifact. Figure [lb indicates that we 
recover a good approximation to the correct, nearly con- 
stant star formation rate in molecular gas at surface den- 
sities from roughly 5 — 100 Mq pc^^. Combined, these 
two effects produce a star formation rate that increases 
superlinearly with total gas content below the H i satu- 
ration threshold and only linearly above it (Figure [TJ;) . 
Third, we recover the return to a superlinear increase of 
star formation rate with total gas content above ^ 100 
Mq pc^^, produced by the increase in molecular cloud 
density in high- pressure environments (Figure [2]) • 

It is worth noting that the observations are subject 
to significant systematic errors in both the gas surface 
density and in the star formation rate, and that this is 
likely to limit the extent of agreement between the data 
and any theoretical model, including this one. The main 
uncertainty in the gas surface density is in the factor 
Xqo used to convert an observed CO intensity into a 
molecular cloud mass. The standard assumption of a 
constant Xco = 2.0 x lO^o cm'^ (K km s~^y^ is likely 
to be accurate to a factor of 2 in spiral galaxy envi- 
ronments, where its value has been calibrated against 



other methods of estimating mass (|Blitz et al.l[20nl . 
In low metallicity dwarf galaxies, however, there is ev- 
idence that the assumption of a fixed Xco may under- 
estimate the tru e molecu lar m ass by as much as a fac- 
tor of ~ 10 (jLerov et alll2007l ). while in starburst sys- 
tems l^jnaY_ovBrestim^e_tl^ mass by factors of up to 
~ 5 (jPownes fc Solomonlll998D . The variation of Xco in 
starbursts has been at least approximately accounted for 
in the high Sg systems shown in Figure [21 but the uncer- 
tainty there is probably larger than for spiral galaxies. 
For dwarf systems at low metallicity, on the other hand, 
the values of Sg and in Figures [1] and [2| have for 
the most part been derived assuming a constant XcOj 
and thus the possible systematic underestimate of Sg 
and h^-s not been included. This would tend to shift 
points to the right in Figure [TJ; and Figure 2. The ef- 
fect in Figure [TJa will be minimal, because most of the 
galaxies for which this effect is significant fall below the 
completeness limit in any event. 

Uncertainties in the star formation rates come from 
a combination of uncertainties in dust corrections and 
in the stellar initial mass function (IMF). Comparing 
star formation rates in the THINGS sample based on 
FUV plus 24 fim emission to those based on Ha, or 
Ha plus 24 /im emission, suggests uncertainties below 
the factor of ~ 2 level. Comparison to the star for- 
mation rate in the Milky Way inferred from radio cat- 
alogs o f H II regions sugg ests a slightly larger uncer- 
tainty: iMcKee fc Williams! ([l997) infer a star formation 
timescale in the molecular gas of (SFRff/iff)~^ = 300 
Myr from this technique, compared to 2 Gyr for the av- 
erage of the THINGS sample. The origin of this discrep- 
ancy is unclear, but it suggests that significant caution 
is warranted in interpreting the star formation rates in- 
ferred from observations. 

4. DISCUSSION 
4.1. Relation to the Previous Work 

It is important to understand how our results here re- 
late to previous work on star formation laws, both phe- 
nomenological and theoretical. Our new star formation 
law, equation (jlOp . is more complex than the simple 
iKennicuttI (|1998[) powerlaw oc S" with n = 1.4±0.15, 
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Fig. 2. — Star formation rate surface density S» as a function of total gas surface density Eg. Lines and contours are tlie same as in 
Figur e [T] Other points are a compilation of literature data from BO^. We show individual apertures in M51 {black dots, Kennicu tt et al.l 
[20071) azimuthal averages (blue circles) in NGC4736 and NGC5055 l lWone fc Blitz 200^), NGC6946 {Crosthwaite & Turner 200^), and 
M51 dSchuster et al.lp007^ . and global averages for starbursts (open green triangles, Kcimicutt 1998), normal spirals {filled green trian- 
ales. IKe'nnicut^ll998l) . and low surface brightness galaxies {yellow diamonds, ^Wydcr ct al. 2009). The gray arrows and labels indicate 
schematically the dominant physical process responsible for setting the slope in each region. 



and the related observational correlation between the CO 
and infrared luminosities of galaxies (assumed to be prox- 
ies for the molecular gas mass and the sta r formation 
rate, respectiv e ly) re ported by iGreve et al.l ((2005) and 
iRiechers et al.l ()2006l ). Such increased complexity is de- 
manded by improving observations, since the data shown 
in Figures [T] and [2] clearly cannot be fit by single pow- 
erlaws. There are very distinct regions where the rela- 
tionship between and Sg is steep, flattens, and then 
steepens again. That said, our mo del is in fact co nsis- 
tent with the data set on which the iKennicuttI ()1998 ^ law 
is based - these are shown by the green open and filled 
triangles in Figure [2]- and these data in turn fall within 
the contours of the THINGS observations. Thus there is 
no inconsistency in the data themselves. 

Wh y then is it possible to fit the data of the lKennicuttj 
(|1998D sample with a single powerlaw? The answer comes 
partly from the fact that the data are averaged over entire 
galaxies, which introduces significant scatter compared 
to the more recent data that is resolved to sub-kpc scales. 
The primary reason for the single powerlaw Kennicutt fit 
with an index ol n k, 1.4, however, is that most of the 
dynamic range in Eg that gives rise to the index of 1.4 

comes from galaxies with Eg ^ 100 Mq pc~^. The same 
is true of the observed correlation between CO and IR 
luminosities in galaxies, where most of the dynamic range 
in the sample comes from starbursts with large surface 
densities and star formation rates. In this regime our 
predicted law, equation p^ . also reduces to a simple 



powerlaw E* cx Eg '^^; the index 1.33 is well within the 
error bars in Kennicutt's fit. Conversely, a fit to only the 
normal galaxies in Kennicutt's sample produces a much 
steeper best-fit index of n = 2.47 ± 0.39, consistent with 
the steeper slope we predict for normal galaxies due to 
the dependence of the H2 fraction on Eg. Thus the classic 
single powerlaw star formation law is in part an artifact 
of fitting a single powerlaw betw een n ormal galaxies and 
star bursts, a poin t also made bv lGao fc Solomon (,2004| ) 
and lWvder etall ((2009h . 

On the theoretical side, the classical explanation for 
the observed E, cx Eg and Lir, oc L^q correlations with 
n K, m « 1.5 is that the star formation rate should 
be pr oportional to the gas mass divid e d by the free-fall 
time (iMadord 119771 ; lElmegreenI 11994 IKennicuttI Il998t 
lElmegreenl 12002). The latter varies as density to the 
— 1/2 power, so in a parcel of gas of density p the star 
formation rate per unit volume varies as p^^"^ . The model 
we propose here is entirely consistent with this basic pic- 
ture, and with earlier explanations for correlations be- 
tween star f ormatio n r ates and g as masses and surface 
densities bv IKM05I and iKrumholz fc Thompson (2007.), 
despite the fact that we obtain a nearly linear correlation 

between E, and Eh^ for Eg < 100 Mq pc "2. 

The reason we can obtain a linear correlation be- 
tween E, and Eua is that, in any theory in which the 
star formation timescale is proportional to the free-fall 
time, the star formation rate per unit area will vary as 
E* oc Efja/^ff oc EjjaP'^^^- Similarly that the total star 
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formation rate in a galaxy will vary as cx MnaP^^^, 
where Mn^ is the total mass of molecular gas. The in- 
dices n and m for the areal and total star formation laws 
therefore depend on how p varies with and A/jja, 
respectively. The most common assum ption in previous 
work has been to adopt p cx Eg (e.g. iElmegreenI l2002t ) 
and p (X AfHa , but this assumption is independent of the 
basic idea that the star formation rate varies inversely 
with the free-fall time. Indeed, if the characteristic den- 
sities of star-forming clouds do not vary from galaxy to 
galaxy, as is approximately (though not exactly) the case 
in our model for galaxies with Eg < 85 A/©, then we 

expect a linear relationship between E* and Eua , even 
though the underlying volumetric s tar formation law is 
nonlinear. iKrumholz fc Thompson! (2007) obtain a lin- 
ear relationship between the luminosities of galaxies in 
the infrared and in HCN(1 0) emission for much the 
same reason. They point out that the high critical den- 
sity of the HCN line ensures that it traces gas at similar 
densities in every galaxy, at least up to the most extreme 
starbursts. Since p is constant, the star formation law is 
linear. In our model here p is roughly constant in galax- 
ies with Eg < 85 Mq pc~^ for a different reason - GMC 
densities are unaffected by galactic environment because 
GMC pressures greatly exceed mean ISM pressures - but 
the end result, a linear star formation law, is the same. 

3 /2 3/2 

Conversely, a scaling law E, oc Ejj^ or Ljr oc L^q is 
expected when the characteristic density of the molecu- 
lar gas scales close to linearly with its surface density or 
total mass. In our model such a scaling appears in the 
starburst regime because ISM pressure becomes impor- 
tant in setting GMC internal densities. 

4.2. Predictions for Future Observations 

Our model makes distinct predictions that can be 
tested against future observations. One obvious predic- 
tion is the upward kink in the relationship between E* 
and Eh2, or E* and Eg, seen at Ehs « Eg « 100 Mq 
pc^^ (Figures la, lb, and 2). This kink in our model is 
caused by the transition between molecular cloud surface 
densities being determined by internal regulation and 
being determined by external pressure. A spatially re- 
solved survey such as THINGS targeting the nearest cir- 
cumnuclear starburst galaxies such as M64, which reach 
surface densities in the r ange ~ 50 — 1000 Mq pc~^ 
(jRosolowskv fc Blit3l2005[ ) , should reveal this kink quite 
clearly, perhaps even within a single galaxy. The one 
caveat is that locating the kink depends on being able 
to measure molecular surface densities accurately, which 
in turn depends on our knowledge of the X factor used 
to convert CO intensity to surface density. Uncertainties 
in its value translate directly into uncertainties in the 
location of the kink in the star formation law. 

A second testable prediction comes from our predicted 
metallicity-dependence of the relationship between E* 
and Eg (Figure Ic). Our model curves span the range 
of metallicities covered by the THINGS sample, and, 
as of the date of this paper's submission, a version of 
the data shown in Figure 1 binned by metallicity was 
not available. However, such metallicity binning should 
reveal that low metallicity galaxies have systematically 
lower star formation rates at fixed Eg for galaxies when 



Eg ^ 10 Mq pc~^. Because there is degeneracy between 
the metallicity Z' and the clumping factor c, ideally this 
prediction would be tested with higher resolution obser- 
vations that reduce or remove the need to adopt a clump- 
ing factor to account for the structure of the gas on scales 
unresolved by the observations. Surveying a large sam- 
ple of entire galactic disks with ~ 100 pc resolution, so 
that c « 1, is probably prohibitively expensive with the 
current generation of telescopes. However, even a ran- 
dom sample of sightlines through galaxies at varying Eg 
and metallicity, with ~ 100 pc resolution, should be suf- 
ficient to test for a metallcity-dependent correlation be- 
tween star formation rate and gas surface density such 
as the one we predict. 

5. SUMMARY 

We have shown that the observed relationship between 
the star formation rate and the atomic and molecular 
content of galaxies can be explained by a simple model, 
whose elements are summarized by the regions labelled 
in Figure [21 First, self-shielding of hydrogen determines 
the amount of gas in molecular form. This imposes a 
characteristic gas surface density of ^ 10 /cZ' Mq pc~^ 
for the transition from atomic to molecular, where c is 
the factor by which the gas surface density is increased 
due to clumping unresolved by the observations and Z' is 
the metallicity relative to solar. Second, once molecules 
do form, molecular clouds reach a surface density of 
roughly 85 M© pc~^ independent of galactic environ- 
ment. This behavior can be understood as arising from 
the fact that molecular clouds are overpressured relative 
to their surroundings, so they must be regul ated by in- 
ternal processes, mos t likely H ii regions ((Matzneij|2002t 
IKrumholz et al.ll2006[ ). that do not depend on metallicity 
or other large-scale galaxy properties. The constant sur- 
face density imposes a roughly constant volume density 
and free-fall time on all molecular gas. The exception 
to this is galaxies where the mean galactic surface den- 
sity is ^ 100 Mq pc^^, in which the ambient pressure 
is high enough to force GMC densities to rise along with 
galactic surface density in order to keep the clouds in 
pressure balance. Third, once formed molecular clouds 
convert themselves into stars at a nearly universal rate 
of ~ 1% of the mass per free-fall time as a result of tur- 
bulent regulation. Together these effects produce a total 
gas star formation law that is superlinear at low galactic 
column density (due to increasing molecular fraction), 
linear or slightly sub-linear at intermediate column (due 
to the invariance of molecular cloud surface densities and 
the weak dependence of GMC masses on galactic prop- 
erties), and superlinear again at high column (due to the 
breakdown of this invariance at high galactic pressures). 

It is worth noting that our model does not make any 
explicit reference to galactic-scale processes such as spi- 
ral shocks, gravitational instability, supernova feedback, 
or cloud-cloud collisions. In a sense all of this physics 
is "upstream" of our theory: processes such as these are 
almost certainly responsible for determining the distri- 
bution of gas surface density within a galaxy. Our model 
addresses the next step of how, once large-scale processes 
assemble the gas, some fraction of it forms molecular 
clouds and then turns into stars. We have therefore sep- 
arated the problem of star formation into two parts, and 
provided a tentative solution for one of them: galactic- 
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scale processes determine Sg, but the physics responsi- 
ble for determining the star formation rate thereafter is 
purely local, and can be understood without reference to 
galactic-scale behavior. Our model shows that a much 
of the recent observational work on star formation can 
be understood in terms of a simple model for that local 
process. 
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